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Thermally activated flux flow and vortex glass transition of recently discovered SmFeAsOo.gFo.i 
superconductor are studied in magnetic fields up to 9.0 T. The thermally activated energy is analyzed 
in two analytic methods, of which one is conventional and generally used, while the other is closer to 
the theoretical description. The thermally activated energy values determined from both methods 
are discussed and compared. In addition, several critical magnetic fields determined from resistivity 
measurements are presented and discussed. 
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I. INTRODUCTION 

The recently discovered FeAs-based superconductors 
inspired study, as their superconducting transition tem- 
peratures and upper critical magnetic fields reach the 
values which are higher than those of MgB2 and com- 
parable to those of the cuprate based superconductors 
(CBS) [J S i i, 1, 1 0, Hi, M EH ■ One of interesting 
characteristics is that they show layered structure with 
conducting layers of FeAs and charge reservoir layers of 
ReO, where Re is a rare earth element Q. This layered 
structure is very similar to that of CBS, and suggests that 
the superconducting behaviors may have similarities to 
those of CBS. 

The vortex dynamics of CBS have been widely stud 
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311 ] . According to the theory, the thermally ac- 
tivated flux flow (TAFF) resistivity is expressed as 
LBJJ) exp(- J c oBVL/T) sinh( JBVL/T) 0, 
22j, where vq is an attempt frequency for 
a flux bundle hopping, L the hopping distance, B the 
magnetic induction, J the applied current density, J c q 
the critical current density in the absence of flux creep, 
V the bundle volume, and T the temperature. If J is 
small enough and JBVL/T <C 1, we have 
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p=(2p c U/T)eM-U/T), 



(1) 



where U = J c qBVL is the thermally activated energy 
(TAE), and p c = vqLB / 'J c q. Equation |T]) simply means 
that the prefactor 2p c lI/T is temperature and magnetic 
field dependent. 

Generally, the TAE of CBS is analyzed by equation ([I]) 
using an assumption that the prefactor 2p c U /T is tem- 
perature independent, and In p(H, T) linearly depends on 
l/T with the form bxp(H,T) = In p (H) - U (H)/T, 
where H is the magnetic field strength, hipo(H) = 
lnpof + U~o(H)/T c [note that In po (H) is the \np(H) value 
for l/T 0], p f the constant, U the TAE for T -> 0, 
and T c the superconducting transition temperature. The 
importance is that the analysis leads to U = Uq(1 — t), 
and the apparent activate energy —d\a pjdT -1 = Uq, 
where t = T/T c [H, H, EE M M, WM M HI . By 



drawing resistivity data in the so-called Arrhenius plot 
with a relation hip(H, l/T), one can easily determine 
Uq(H) with its corresponding slope in a low resistivity 
range. However , U = Uo(l — t) may not be true in real- 
ity, when a local slope of In p vs l/T in the Arrhenius plot 
shows a round curvature. As a result, the corresponding 
apparent activated energy —dh\p(H,T)/dT~ 1 shows a 
sharp increase with decreasing temperature. The phe- 
nomena means U ^ Uq(1 — t), pof ^ const, and the de- 
termination of Uq(H) from the slopes is a problem. In the 
early stage of the discovery of CBS, the experimental ob- 
servation of the abnormal phenomena had been reported 
by Palstra et al [22J without solution. 

Zhang et al [23j suggested that the temperature de- 
pendent of the prefactor in equation ([1]) must be taken 
into account in the analysis. By using this suggestion to 
equation |T]) with U = Uq{1 — t) q , the derivative, 



dlnp/dT- 1 = (1 - T/U)(U - TdU/dT) 

= [U (l-ty-T}[l + qt/(l-t)}, (2) 



where q has a value in the range from 0.5 to 2, and t = 
T/T c . 

In this paper, TAFF resistivity of SmFeAsOo.gFo.i 
(SFAOF) is studied with magnetic fields up to 9.0 T. 
By using the assumption pof = 2p c U/T = const, and 
U = Uo(l—t), the TAFF behaviors were first analyzed in 
the Arrhenius plot. After that, by using the assumption 
that the prefactor 2p c ll /T is temperature dependent, and 
p c is not temperature dependent, the TAFF resistivity is 
analyzed in the other method in which equation ^ is em- 
ployed to determine Uo(H). The Uo{H) determinations 
from both methods are discussed and compared. We sug- 
gest that the second method shall be instead of the first 
in the analysis of TAFF characteristics of other super- 
conductors. In addition, the vortex glass transition, zero 
resistivity temperature, and the temperature dependents 
of the critical fields determined from different resistivity 
criteria values in the superconducting transition regime 
are presented. 



2 



II. EXPERIMENTS 

The SFAOF sample was prepared by a high-pressure 
synthesis method. SmAs powder (pre-sintered) and 
As, Fe, Fe2C>3, FeF2 powders (the purities of all start- 
ing chemicals are better than 99.99%) were mixed 
together with the nominal stoichiometric ratio of 
SmfOi-ajF^jFeAs, then ground thoroughly and pressed 
into small pellets. The pellet was sealed in boron nitride 
crucibles and sintered in a high pressure synthesis ap- 
paratus under the pressure of 6.0 GPa and temperature 
of 1250° C for two hours. The x-ray diffraction analysis 
showed that the main phase is LaOFeAs structure with 
some impurity phases [2], |3J|. It was cut into a rectan- 
gular shape with dimensions of 4.20 mm(length)x 1.60 
mm(width) x 1.08 mm(thickness). The standard four 
probe technique was used for resistivity p(T, H) mea- 
surements. Bipolar pulsed dc current with an amplitude 
of 5.0 mA (corresponding to the current density about 
^0.29 A/cm 2 ) was applied to it. The measurements were 
performed on a physical property measurement system 
(PPMS, Quantum Design) with the magnetic field up to 
9 T. From the zero field p(T, 0) data, we find that the 
superconducting transition width is about 1.7 K (defined 
by the superconducting transition from 10% to 90% of 
the normal resistivity), and the zero resistant temper- 
ature T c q is 52.2 K (determined by the criterion of 0.1 
/ifi-cm). 



III. RESULTS AND DISCUSSION 

Figure Dla)-(c), respectively, show p(T,H), 
-d\np(T,H)/dT-\ and dT/ In p(T, H) data with 
different symbols for 0.0, 0.5, 1.0, 3.0, 5.0, 7.0, and 
9.0 T fields. One shall notice that the suppression 
of the superconductivity at 9.0 T field is comparable 
to that of low anisotropic CBS, e.g. YBa2Cu307_a 
(YBCO), and suggests that the anisotropy is similar 
to YBCO. The so-called apparent activated energy 
—d In p(T, HydT -1 data were generally used to analyze 
TAE (23, Hi, and the derivative dT/ In p{T,H) data 
is always used to analysis the vortex glass transition 
boundary [H, HH, [H, [35|. Below, we will present our 
analysis. 

Figure [5] shows the Arrhenius plot of its resistivity in 
magnetic fields of poH = 0.0, 0.5, 1.0, 3.0, 5.0, 7.0, and 
9.0 T with different symbols. The solid lines plot lin- 
ear regressions in the low resistivity range. According 
to the conventional analysis method, these linear fits are 
based on the assumption U(T,H) — Uq(H)(1 — t), where 
each Uo(H) was determined by using each slope. Note 
that there is a cross point where all the linear fitting 
lines approximately focus together except for the line for 
zero field, and this point leads to T cr w 53.8 K. The in- 
set shows \npo(Uo) data and the corresponding solid line 
with the linear fitting In po = hipof+Uo/T c . From the fit- 
ting, we determined pof = 6.668 mi7-cm, and T c « 53.8 K 




FIG. 1: (a)-(c), respectively, show p(T,H), 

~dln P (T,H)/dT-\ and dT/dlnp(T,H) data in mag- 
netic fields of p H = 0.0, 0.5, 1.0, 3.0, 5.0, 7.0, and 9.0 
T with different symbols. The solid lines in (a) and short 
horizontal lines in (b) correspond to the regressions with the 
fitting parameters of Uo(H) and pof determined from the 
first analytic method. The dashed lines are regression curves 
with the fitting parameters Uo(H) and p c (H) determinded 
from the second analytic method (see text). The dashed lines 
in (c) plot linear fittings. 




3x10 2 2.0x10 2 2.2x10 2 2.4x10 2 2.6x10 2 
1/F(1/K) 

FIG. 2: Arrhenius plot of p(T,H) for p H = 0.0, 0.5, 1.0, 3.0, 
5.0, 7.0, and 9.0 T. The solid lines are plot linear fitting for 
low resistivity. The dashed lines are regressive curves using 
the second analytic method. The inset shows \n P o(Uo) data 
determined by using the first analytic method, where the solid 
line is a linear fit for the data. 
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which coincides with the value of T cr . Using the pof, T c , 
and Uq(H) data, we regressed the p(T, H) data as shown 
in figure QJa) with solid lines. Since the assumption of 
U(T,H) = U (H)(l-t) leads to -dlnp^H)/^- 1 = 
Uq, we present Uq{H) data in figure [TJb) with horizon- 
tal solid lines where each of them has a limited length. 
Each length covers the temperature interval which corre- 
sponds to the interval of the reciprocal temperature for 
determining Uq(H) in the Arrhenius plot. For decreasing 
temperature, note that each — d\np(T, J7)/9T _1 dataset 
approximately intersects its horizontal Uq(H) line cen- 
ter with a divergent trend in the temperature interval. 
This means that each Uq{H) approximates to the av- 
erage value of its —d\np(T,H)/dT~ 1 in the tempera- 
ture interval. The similar divergent trends were also ob- 
served in YBCO [22L |23(| . The analysis indicates that the 
TAE determined from the conventional method may have 
problems in characteristics, and the method suggested by 
Zhang et al [23| ought to be taken into consideration. 

In figure D^b), one may have noticed that each 
— d\np(T, H")/dT~ l set can be divided into five regimes: 
(i) the normal state regime at high temperature, 
where —dlrLp(T,H)/dT~ 1 is almost temperature-and 
magnetic- field-independent; (ii) the primary supercon- 
ducting transition regime, where — dh\p(T, H)/dT~ 1 
quickly increases and the resistivity begins sharply de- 
creasing (see figureQJa)); (iii) the platform regime, where 
— d\np(T, _ff)/<9T -1 for each field measurement show a 
step structure (see figure D^b)); (iv) the second sharply 
increasing regime, where —dhip(T,H)/dT~ 1 quickly in- 
creases into a high value range and dT/d\np data show 
a linear type (see figure HJc)); (v) the strong fluctua- 
tion regime, where —d\np(T,H)/dT~ 1 curve shows an 
irregular shape due to the resistivity reaching the low 
measurable range. 

For analyzing TAFF behaviors with the second ana- 
lytic method, the first thing is to validate one or two 
regimes which relate to the TAFF behaviors. Appar- 
ently, the data in regimes (i) and (ii) do not relate to 
TAFF behaviors, as the data are in the normal state 
in regime (i) and in the flux flow regime in regime (ii). 
The data show a platform structure and possibly relate 
to TAFF behaviors in regime (iii). Note that resistivity 
data in the regime are only about one order of magni- 
tude less than that in regime (i) (where the resistivity is 
in the normal state) [see figure DJa) an d (b)]. Therefore, 
we conclude that the data in regime (iii) are not in the 
TAFF regime [22J . In regime (iv) , the resistivity is about 
2 and 3 orders of magnitude less than that in regime (i). 
The resistivity in this range, as suggested by Palstra et 
al [22j], relates to TAFF behavior. According to the vor- 
tex glass transition theory, the linear (ohmic) resistivity 

where T g 
Ac- 



ought to linearly vanish in a form p oc (T — T„ 
is the vortex glass transition temperature [32 
cordingly, dT/d\np oc (T - T g ) 0, H|. In figure QJb), 
one will easily find that dT/ din p(T, H) curves show ap- 
proximately linear curvatures in regime (iv). The data in 
the regime (v) show irregular trend, as the resistivity is 




foH (T) 

FIG. 3: (a) The circles show Uo(H) extracted from the slopes 
in Arrhenius plot in figure [2] The downward triangles show 
Uo(H) data determined by using the relation —d\n p/dT^ 1 = 
[U (l-t) q -T]{l+qt/(l-t)], where q = 2. The stars represent 
corresponding p c data. 



going to a deeply superconducting state where the resis- 
tivity may be dominated by non-ohmic characteristics; 
besides, the measuring voltmeter reached its low mea- 
suring limitation in experiments. The analysis concludes 
that the TAFF resistivity is in regime (iv). 

Figure [3] shows Uq{H) determined by the first ana- 
lytic method with circles and by the second with trian- 
gles. The stars show corresponding p c (H) for the second 
method. All the dash lines in figures [T^a) and (b), and 
figure^ are regression curves with the regression param- 
eters of Uo(H) and p c {H) in figure [3l In the analysis of 
these data with the second method, we first derived U$ 
with equation as Uq is the only free parameter in the 
equation except for q. We found that the energy relation, 
U{T,H) = U Q (1 - t) q with q = 2, leads to good consis- 
tency to the experimental data, where t = T/T c , and 
T c = 53.8 K. After determination of each Uq(H), each 
p c (H) can be easily determined by fitting equation (jTJ) . 
One will find that all the regressions (dashed lines) in fig- 
ure da) , figure []Jb) , and figure [2] are in good agreement 
with experimental data and confirm the correctness of 
the analytic method. Note that, although each tempera- 
ture interval in regime (iv), which we used to determine 
Uq (H) with the second analytic method, is somewhat less 
than that we used to determine Uo(H) with the first an- 
alytic method (see figure QJb) and (c)), the regressions of 
the second method are still given better fitting results. 

For magnetic field above 1.0 T, we find that U oc 
_ff~°- 57 (l — t) for the data derived from Arrhenius plot 
with the first method, and U oc H~°"(l — t) 2 with the 
second method. Note that the TAE determined by the 
second method is about one order larger than that deter- 
mined from the first as shown in the figure. The reason 
of the large difference between the two analytic methods 
is due to that the first method not taking into account 
the temperature dependent relation of the prefactor in 
equation (jTJ) , while the second method takes into consid- 
eration. Apparently, the second method is closer to the 
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FIG. 4: Critical magnetic fields. The circles show the H(T g ) 
data, the downward triangles the H(T c o), the squares the 
H(T 10 %), the diamonds the H(T 50 %), and the upward tri- 
angles the H(T 9Q% ). 



theoretical description. 

Figure 2] shows the vortex glass transition boundary 
H(T g ), H(T c0 ), H(T 10% ), H(T 50% ), and H(T 90% ), where 
T g is the glass transition temperature, T c q the zero 
resistance temperature, T 10 % the 10% of the normal 
state resistivity temperature, T 50 % the 50% and T go % 
the 90%. Here, H(T g ) was determined by using the 
dT/dlnp oc (T — T g ) relation [dash lines in figure 
[He)]. We find that the data can be fitted by a rela- 
tion H g = H g0 {\ - Tg/52.05) 2 , where fi H g0 « 201 T. 
The poHgo value may not be true in zero temperature, 
but it leads to a conclusion that the vortex melting field 
is very large for SFAOF superconductor when the tem- 
perature approaches to zero. A similar characteristic can 
be also observed for H(T c0 ) data, since each T c0 (H) is 
slightly higher than that T g {H). The normalized transi- 
tion widths (defined as (1 — Ti %/T 90 %) are comparable 
to YBCO, and thus suggest that the anisotropy of the 
SFAOF is similar to YBCO. The H{T gQ% ) curve shows a 
steep increase indicating that the upper critical magnetic 
field of the superconductor is very high in low tempera- 
ture. One may notice that the transition widths between 
T c a(H) and T W %{H) are comparable to that between 
T 10 %(H) and T 90 %(H), suggesting that superconductor 
may have better analysis results when further improving 



the quality of the superconductor. 

In comparisons, a sample that was prepared in the 
same sample batch was used in ac susceptibility measure- 
ments with magnetic fields up to 7.0 T, and frequencies 
up to 1.11 kHz. We found that T g (H) is somewhat larger 
than the corresponding temperature of the peak in the 
imaginary component in ac susceptibility measurement. 
However, the similar H increasing trends were also found 
in the ac measurement. Besides, we found that the onset 
temperature of the superconducting transition T onset (H) 
of ac susceptibility measurement is somewhat higher than 
that ofT c o(H). Detailed analysis of the ac measurement 
is just ongoing. 



IV. CONCLUSION 

In summary, we analyze the resistive TAFF behavior 
of SFAOF with two theoretical analysis. For the first 
method, pof = const and U = Uo(l — t) were assumed, 
and thus Arrhcnius plot was employed in analysis. This 
method is simple and easy in analysis, but the anal- 
ysis results remain inconsistency to experimental data. 
The second method assumes that the prefactor 2p c lI/T 
is temperature dependent, while p c is not temperature 
dependent. By using the second method, equation ([2]) is 
obtained. The second method results in the regressions 
are in good agreement with experimental data. The TAE 
analysis shows U ~ iJ _1 (l — t) 2 for magnetic field above 
1.0 T. The second method is also simple and easy, since 
only two free parameters in analysis. We suggest that the 
second method shall be instead of the first for the anal- 
ysis of TAFF behaviors of other superconductors. The 
study shows that critical magnetic fields of SFAOF may 
have large values in low temperature and comparable to 
CBS. 



V. ACKNOWLEDGEMENTS 

We are grateful to Prof. H. H. Wen and Dr. H. Yang 
for resistivity measurements and helpful discussion. This 
work has been financially supported by the National Nat- 
ural Science Foundation of China (Grant No. 10874221). 



[1] Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, 

J. Am. Chem. Soc. 130, 3296 (2008). 
[2] Z. A. Ren, W. Lu, J. Yang, W. Yi, X. L. Shen, Z. C. Li, 

G. C. Che, X. L. Dong, L. L. Sun, F. Zhou, and Z. X. 

Zhao, Chin. Phys. Lett. 25, 2215 (2008). 
[3] Z. A. Ren, G. C. Che, X. L. Dong, J. Yang, W. Lu, W. 

Yi, X. L. Shen, Z. C. Li, L. L. Sun, F. Zhou, Z. X. Zhao, 

Europhysics Letters 83, 17002 (2008). 
[4] H. H. Wen, G. Mu, L. Fang, H. Yang, and X. Zhu, Eu- 

rophys. Lett. 82, 17009 (2008). 



[5] G. F. Chen, Z. Li, D. Wu, G. Li, W. Z. Hu, J. Dong, 
P. Zheng, J. L. Luo, N. L. Wang , Phys. Rev. Lett. 100, 
247002 (2008). 

[6] G. F. Chen, Z. Li, D. Wu, J. Dong, G. Li, W. Z. Hu, P. 
Zheng, J. L. Luo, N. L. Wang, Chin. Phys. Lett. 25, 2235 
(2008). 

[7] Clarina de la Cruz, Q. Huang, J. W. Lynn, Jiying Li, 
W. Ratcliff II, J. L. Zarestky, H. A. Mook, G. F. Chen, 
J. L. Luo, N. L. Wang, Pengcheng Dai, Nature 453, 899 
(2008). 



5 



[8] J. Dong, H. J. Zhang, G. Xu, Z. Li, G. Li, W. Z. Hu, D. 

Wu, G. F. Chen, X. Dai, J. L. Luo, Z. Fang, N. L. Wang, 

Europhysics Letters, 83, 27006 (2008). 
[9] F. Hunte, J. Jaroszynski, A. Gurevich, D. C. Larbalestier, 

R. Jin, A. S. Sefat, M. A. McGuire, B. C. Sales, D. K. 

Christen, and D. Mandrus, Nature (London) 453, 903 

(2008). 

[10] C. Senatore, M. Cantoni, G. Wu R.H. Liu, X.H. Chen, 

and R. Fhikiger, larXiv:0805.2389l 
[11] Y. Jia, P. Cheng, L. Fang, H. Q. L uo, H. Yang, C. R en, 

L. Shan, C. Z. Gu, and H. H. Wen, larXiv:0806.0532l 
[12] Y. Yeshurun and A. P. Malozemoff, Phys. Rev. Lett. 60, 

2202 (1988). 

[13] A. P. Malozemoff, L. Krusin-Elbaum, D. C. Cronemeyer, 
Y. Yeshurun and F. Holtzberg, Phys. Rev. B 38, 6490 

(1988) . 

[14] C. W. Hagen and R. Griessen, Phys. Rev. Lett. 62, 2857 

(1989) . 

[15] V. B. Geshkenbein, M. V. Feigel'man, A. I. Larkin, and 
V. M. Vinokur, Physica C 162, 239 (1989). 

[16] P. H. Kes, J. Aarts, .J van den Berg, C. J. van der 
Beek and J. A. Mydosh, Supercond. Sci. Technol. 1, 242 
(1989). 

[17] V. M. Vinokur, M. V. Feigel'man, V. B. Geshkenbein, 

and A. I. Larkin, Phys. Rev. Lett. 65, 259 (1990); J. 

Kierfeld, H. Nordborg, and V. M. Vinokur, Phys. Rev. 

Lett. 85, 4948 (2000). 
[18] G. Blatter, M. V. Feigel'man, V. B. Geshkenbein, A. I. 

Larkin, and V. M. Vinokur, Rev. Mod. Phys. 66, 1125 

(1994), (and the references therein). 
[19] E. H. Brandt, Rep. Prog. Phys. 58, 1465 (1995), (and 

the references therein). 
[20] D. Shi, High- Temperature Superconducting Materials 

Science and Engineering New Concepts and Technology, 

(Pergamon, New York, 1995). 
[21] M. Tinkham, Introduct. to Superconduct., McGraw-Hill, 



New York, 1996, (and the references therein). 
[22] T. T. M. Palstra, B. Batlogg, L. F. Schneemeyer and J. 

V. Waszczak, Phys. Rev. Lett. 61, 1662 (1988); T. T. M. 

Palstra, B. Batlogg, R. B. van Dover, I. F. Schneemeyer, 

and J. V. Waszczak, Phys. Rev. B 41, 6621 (1990). 
[23] Y. Z. Zhang, H. H. Wen, Z. Wang, Phys. Rev. B 74, 

144521 (2006). 

[24] Y. Z. Zhang, Z. Wang, X. F. Lu, H. H. Wen, J. F. de 
Marneffe, R. Deltour, A. G. M. Jansen, and P. Wyder, 
Phys. Rev. B 71, 052502 (2005). 

[25] J. Figueras, T. Puig, and X. Obradors Phys. Rev. B 67, 
014503 (2003) 

[26] Cao Xiaowen, Wang Zhihe, and Xu Xiaojun Phys. Rev. B 
65, 064521 (2002); Xu Xiaojun, Fu Lan, Wang Liangbin, 
Zhang Yuheng, Fang Jun, Cao Xiaowen, Li Kebin, and 
Sekine Hisashi Phys. Rev. B 59, 608 (1999). 

[27] Y. Z. Zhang, R. Deltour, and Z. X. Zhao, Phys. Rev. 
Lett. 85, 3492 (2000). 

[28] M. Andersson, A. Rydh, and 6. Rapp Phys. Rev. B 63, 
184511 (2001). 

[29] D. Ravelosona, J. P. Contour, and N. Bontemps Phys. 

Rev. B 61, 7044 (2000). 
[30] S. N. Gordeev, A. P. Rassau, R. M. Langan, P. A. J. de 

Groot, V. B. Geshkenbein, R. Gagnon, and L. Taillefer 

Phys. Rev. B 60, 10477 (1999). 
[31] H. C. Yang, L. M. Wang, and H. E. Horng Phys. Rev. B 

59, 8956 (1999). 
[32] M.P.A. Fisher, Phys. Rev. Lett. 62, 1415 (1989). 
[33] D.S. Fisher, M.P.A. Fisher, and D.A. Huse, Phys. Rev. 

B 43, 130 (1991). 
[34] P. Wagner, U. Frey, F. Hillmer, and H. Adrian, Phys. 

Rev. B 51, 1206 (1995). 
[35] Y. Z. Zhang, R. Deltour, J.-F. de Marneffe, Y. L. Qin, 

L. Li, Z. X. Zhao, A. G. M. Jansen, and P. Wyder Phys. 

Rev. B 61, 8675 (2000). 



